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Abstract 

The ATLAS experiment at the LHC has measured the Higgs boson couplings and mass, and 
searched for invisible Higgs boson decays, using multiple production and decay channels 
with up to 4.7 fb”^ of pp collision data at yfs = 7 TeV and 20.3 fb”^ at aA = 8 TeV. In 
the current study, the measured production and decay rates of the observed Higgs boson in 
the yy, ZZ, WW, Zy, bb, rr, and pp decay channels, along with results from the associated 
production of a Higgs boson with a top-quark pair, are used to probe the scaling of the coup¬ 
lings with mass. Limits are set on parameters in extensions of the Standard Model including 
a composite Higgs boson, an additional electroweak singlet, and two-Higgs-doublet models. 
Together with the measured mass of the scalar Higgs boson in the yy and ZZ decay modes, a 
lower limit is set on the pseudoscalar Higgs boson mass of > 370 GeV in the “hMSSM” 
simplified Minimal Supersymmetric Standard Model. Results from direct searches for heavy 
Higgs bosons are also interpreted in the hMSSM. Direct searches for invisible Higgs boson 
decays in the vector-boson fusion and associated production of a Higgs boson with W/Z 
(Z ££, WIZ jj) modes are statistically combined to set an upper limit on the Higgs bo¬ 
son invisible branching ratio of 0.25. The use of the measured visible decay rates in a more 
general coupling fit improves the upper limit to 0.23, constraining a Higgs portal model of 
dark matter. 
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1 Introduction 


The ATLAS and CMS Collaborations at the Large Hadron Collider (LHC) announced the discovery of 
a particle consistent with a Higgs boson in 2012 [1, 2]. Since then, the collaborations have together 
measured the mass of the particle to be about 125 GeV [3-5]. Studies of its spin and parity in bosonic 
decays have found it to be compatible with a 7^ = O'*" state [6-8] . Combined coupling fits of the measured 
Higgs boson production and decay rates within the framework of the Standard Model (SM) have found 
no significant deviation from the SM expectations [4, 9, 10]. These results strongly suggest that the newly 
discovered particle is indeed a Higgs boson and that a non-zero vacuum expectation value of a Higgs 
doublet is responsible for electroweak (EW) symmetry breaking [11-13]. The observed CP-even Higgs 
boson is denoted as h throughout this paper. 

A crucial question is whether there is only one Higgs doublet, as postulated by the SM, or whether the 
Higgs sector is more complex, for example with a second doublet leading to more than one Higgs boson of 
which one has properties similar to those of the SM Higgs boson, as predicted in many theories beyond the 
Standard Model (BSM).* The “hierarchy problem” regarding the naturalness of the Higgs boson mass, the 
nature of dark matter, and other open questions that the SM is not able to answer also motivate the search 
for additional new particles and interactions. Astrophysical observations provide strong evidence of dark 
matter that could be explained by the existence of weakly interacting massive particles (see Ref. [14] 
and the references therein). If such decays are kinematically allowed, the observed Higgs boson [1, 2] 
might decay to dark matter or other stable or long-lived particles which do not interact significantly with a 
detector [15-20]. Such Higgs boson decays are termed “invisible” and can be inferred indirectly through 
final sfafes wifh large missing fransverse momenfum. The Higgs boson may also decay fo particles fhaf do 
inferacf significanfly wifh a defector, such as gluons fhaf produce jefs, resulting in final sfafes fhaf cannof 
be resolved due fo fhe very large backgrounds. These decays and final sfafes are fermed “undefecfable”. 

This paper presenfs searches for deviafions from fhe rafes of Higgs boson production and decay predicfed 
by fhe SM, including bofh fhe visible and invisible decay channels, using ATLAS dafa. Simulfaneous 
fifs of multiple production and decay channels are performed affer fhe removal of overlaps in fhe even! 
selection of differenf analyses, and correlations befween fhe sysfemafic uncerfainfies are accounted for. 
The dafa are interpreted in various benchmark models beyond fhe SM, providing indirecf limifs on fhe 
BSM paramefers. The limifs make differenf assumpfions fhan fhose obfained by direcf searches for heavy 
Higgs bosons and invisible Higgs boson decays. 

An overview of fhe experimenfal inpufs is given in Secfion 2, and fhe analysis procedure is described in 
Section 3. The scaling of fhe couplings wifh mass is probed in Secfion 4. The measuremenfs of visible 
Higgs boson decay rales are used to derive limifs on model paramefers in four represenlafive classes 
of models: Minimal Composife Higgs Models (MCHM) in Secfion 5, an additional elecfroweak singlel 
in Secfion 6, fwo-Higgs-doublef models (2HDMs) in Secfion 7, and fhe “h” Minimal Supersymmelric 
Slandard Model (hMSSM) in Secfion 8. The resulls from direcf searches for heavy Higgs bosons are also 
inferprefed in fhe hMSSM in Secfion 8. The combination of direcf searches for invisible Higgs boson 
decays is discussed in Section 9.1, and fhe combination of all visible and invisible Higgs boson decay 
channels is described in Section 9.2. This is used logefher wifh fhe visible decays to consfrain a Higgs 
porlal model of dark mailer in Secfion 9.3. Finally, Secfion 10 is devoled fo fhe conclusions. 


* The observed CP-even Higgs boson, denoted as h in this paper, is taken to be the lightest Higgs boson, and only heavier 
additional Higgs bosons are considered. 
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2 Experimental inputs 


For the determination of the couplings in the visible Higgs boson decay channels, the experimental in¬ 
puts include search results and measurements of Higgs boson decays: h^yy [21], h^ZZ*^A{ [22], 
h^WW*^Mv [23, 24], h^Zy [25], h ^ bb [26], h ^ tt [27], and h ^ fi/j [28] {£ - Search 
results from tth associated production with h^yy [29], h ^ bb [30], and final states with multiple 
leptons [31] are included. In addition, the constraints on the Higgs boson invisible decay branching ratio 
use direct searches for Higgs boson decays to invisible particles in events with dileptons or dijets with 
large missing transverse momentum, These inputs include the search for a Higgs boson, produced 

through vector-boson fusion (VBF) and thus accompanied by dijets, that decays invisibly and results in 
missing transverse momentum (VBF —> 77 +^™**) [32]; the search for a Higgs boson, which subsequently 
decays invisibly, produced in association with a Z boson that decays to dileptons {Zh -^ £{ + [33]); 

and the search for a Higgs boson, which afterwards decays invisibly, produced together with a IF or Z 
boson that decays hadronically (WIZh 77 - 1 - [34]). These searches are based on up to 4.7 fb“^ of 

pp collision data at = 7 TeV and up to 20.3 fb“^ at = 8 TeV. 

Each measurement or search classifies candidate events into exclusive categories based on the expected 
kinematic properties of different Higgs boson production processes. This both improves the sensitivity 
and enables discrimination between different Higgs boson production modes. Each search channel is 
designed to be mostly sensitive to the product of a Higgs boson production cross section and decay 
branching ratio. The combination of the visible decay search channels is used [10] to determine the 
couplings of the Higgs boson to other SM particles. The input analyses, their results, and small changes 
to them applied for use in this combination are described there. 

Direct searches for additional heavy Higgs bosons (H, A, and H-) are not used in the fits discussed here, 
but their results are interpreted in the hMSSM benchmark model for comparison. 


3 Analysis procedure 


In the benchmark models considered, the couplings of the Higgs boson to fermions and vector bosons 
are modified by funcfions of fhe model paramefers. In all cases, if is assumed fhaf fhe modificalions 
of fhe couplings do nof change fhe Higgs boson producfion or decay kinemafics significanfly. Thus fhe 
expecfed rafe of any given process can be obfained fhrough a simple rescaling of fhe SM couplings and 
no accepfance change due fo kinemafics in each BSM scenario is included. A simulfaneous fil of fhe 
measured rales in mulliple producfion and decay modes is used lo conslrain fhe BSM model paramefers. 
The Higgs boson mass was measured by ATEAS fo be m/, = 125.36 ± 0.37 (sfaf) ±0.18 (sysf) GeV [3]. 
The besl-fil value is used fhroughouf fhis paper; fhe uncerfainfy on fhe mass is nof included. 


The sfafisfical frealmenf of fhe dafa is described in Refs. [35-39]. Confidence inlervals use fhe lesl slalisfic 
ta = -21n A(Qr), which is based on fhe profile likelihood ratio [40]: 


ACrr) 


L(a , 6{a)) 
L{a, 6) 


( 1 ) 


The likelihood in Eq. (1) depends on one or more parameters of inferesf or, such as fhe Higgs boson 
producfion times branching ralio strengfh p, fhe mass m/„ and coupling scale factors /c,. Sysfemafic 
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uncertainties and their correlations [35] are modelled by introducing nuisance parameters 9 centred at 
their nominal values. For the visible decay channels, the treatment of systematic uncertainties is the 
same as that used in Higgs boson coupling measurements [10]. For the invisible decay channels, the 
expected event counts for the signals, backgrounds and control regions are taken from Monte Carlo (MC) 
predictions or data-driven estimations as described in Refs. [32-34]. The nuisance parameters for each 
individual source of uncertainty are applied on the relevant expected rates so that the correlated effects of 
the uncertainties are taken into account. 

The single circumflex in the denominator of Eq. (1) denotes the unconditional maximum-likelihood es¬ 
timate of a parameter. The double circumflex in the numerator denotes the “profiled” value, namely the 
conditional maximum-likelihood estimate for given fixed values of the parameters of interest a. 

For each production mode j and visible decay channel k, jj. is normalised to the SM expectation for that 
channel so that ji = \ corresponds to the SM Higgs boson hypothesis and = 0 to the background-only 
hypothesis: 

Cr;,SM X SM 

where ctj is the production cross section, BR^; is the branching ratio, and the subscript “SM” denotes their 
SM expectations. 

For the invisible decay mode, ju is the production cross section for each production mode j times the 
invisible decay branching ratio BRinv, normalised to the total SM rate for the production mode in question: 

jj. = —^ X BRinv ■ (3) 

C j,SM 

Thus the SM is recovered at // = 0 when BRinv= 0. 

Other parameters of interest characterise each particular scenario studied, including the mass scaling 
parameter e and the “vacuum expectation value” parameter M for the scaling of the couplings with mass 
(Section 4), compositeness scaling parameter ^ for the Higgs boson compositeness models (Section 5), 
squared coupling k'^ of the heavy Higgs boson in the electroweak singlet model (Section 6), cos(yS - a) 
and tan/3 for the 2HDM (Section 7), pseudoscalar Higgs boson mass ma and tanyS for the hMSSM model 
(Section 8), and Higgs boson invisible decay branching ratio BRinv for the studies of Higgs boson invisible 
decays (Section 9). 

The likelihood function for the Higgs boson coupling measurements is built as a product of the likelihoods 
of all measured Higgs boson channels, where for each channel the likelihood is built using sums of signal 
and background probability density functions in the discriminating variables. These discriminants are 
chosen to be the yy and fi/u mass spectra for h^yy [21] and h ^ fiji [28] respectively; the transverse 
mass, mj, distribution^ for VFVF*—> £v£v [23, 24]; the distribution of a boosted decision tree (BDT) 

response for h ^ tt [27] and h ^ bb [26]; the At mass spectrum and a BDT in the h^ZZ*^4£ 
channel [22]; the distribution for the VBF —> jj + E™®* [32], Zh ^ it + E™*® [33], and WIZh 
jj + [34] channels. The distributions are derived primarily from MC simulation for the signal, and 

both the data and simulation contribute to them for the background. 

The couplings are parameterised using scale factors denoted x,, which are defined as the ratios of the 
couplings to their corresponding SM values. The production and decay rates are modified from their SM 

^ The transverse mass mj is defined as: mj = yj(Ej + Pj f- — where Eij = , Pj ip^) is the vector 

sum of the lepton (neutrino) transverse momenta, and iPj) is its modulus. 
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expectations accordingly, as expected at leading order [41]. This procedure is performed for each of the 
models probed in Sections 4-9, using the coupling parameterisation given for each model. For example, 
taking the narrow-width approximation [42, 43], the rate for the process gg ^ h ^ ZZ* At relative to 
the SM prediction can be parameterised [41] as: 

^ ctxBR ^ 

(Tsm X BRsm 


Here Kg is the scale factor for the loop-induced coupling to the gluon through the top and bottom quarks, 
where both the top and bottom couplings are scaled by Kf, and kz is the coupling scale factor for the Z 
boson. The scale factor for the total width of the Higgs boson, k^, is calculated as a squared effective 
coupling scale factor. It is defined as fhe sum of squared coupling scale faclors for all decay modes, 
each weighted by the corresponding SM partial decay width T™ [41]: 



J jj 


rSM 
^ 1. 


(5) 


where T™ is the SM total width and the summation runs over WW, ZZ, yy, Zy, gg, tt, bb, cc, ss, tt, 
and nji. The present experimental sensitivity to Higgs boson decays to charm and strange quarks with 
the current data is very low. Therefore the scale factors of the corresponding couplings are taken to be 
equal to those of the top and bottom quarks, respectively, which have the same quantum numbers. The 
couplings to the first-generation quarks (up and down) and the electron are negligible. 


In most of the models considered (Sections 4-8), it is assumed that no new production or decay modes 
beyond those in the SM are kinematically open. In addition, the production or decays through loops 
are resolved in terms of the contributing particles in the loops, taking non-negligible contributions only 
from SM particles. For example, the W boson provides the dominant contribution to the h ^ yy decay 
(followed by the top quark), such that the effective coupling scale factor Ky is given by: 


Ki{Kb, Kt, Kr, Kw) ^ 


X KiKj ■ T^y 

y f'-^ 


( 6 ) 


where Fyy is the contribution to the diphoton decay width due to a particle loop (/ = j) or due to the 
interference between two particles (/ j), and where the summations run over the W boson, top and 

bottom quarks, and tau lepton. Contributions from other charged particles in the SM are negligible. The 
destructive interference between the W and top loops, as well as the contributions from other charged 
particles in the loops, are thus accounted for. Similarly, for the loop-induced h ^ Zy and gg ^ h 
processes the effective coupling scale factors are given by: 


KZyiXb, Kt, Kr, Kw) 


Z KiKj ■ 

U(i>j) 





(V) 


Kn{Kb,Kt) ^ 


cr” , + kI' cr^^, + KfKb ■ , 

ggh b ggh ‘ " ggh 


cr' 


-I- cr^^, -I- cr‘ 


ggh ggh 


tb 


( 8 ) 
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where ^ are the respective contributions to the gluon fusion cross section from a top loop, 

bottom loop, and the interference of the top and bottom loops. 

In the searches for Higgs boson decays to invisible particles discussed in Section 9, it is assumed that 
there are no new production modes beyond the SM ones; however, the possibility of new decay modes 
is left open. The couplings associated with Higgs boson production and decays through loops are not 
resolved, but rather left as effective couplings. 

Confidence intervals are extracted by taking ta to follow an asymptotic distribution with the corres¬ 
ponding number of degrees of freedom [40]. For the composite Higgs boson (see Section 5), EW singlet 
(Section 6 ), and invisible Higgs boson decays (Section 9), a physical boundary imposes a lower bound on 
the model parameter under study. The confidence intervals reported are based on the profile likelihood 
ratio where parameters are restricted to the allowed region of parameter space, as in the case of the test 
statistic described in Ref. [40]. This restriction of the likelihood ratio to the allowed region of parameter 
space is similar to the Feldman-Cousins technique [44] and provides protection against artificial exclu¬ 
sions due to fluctuations into the unphysical regime. However, the confidence interval is defined by the 
standardT^ cutoff, leading to overcoverage near the physical boundaries as demonstrated by toy examples. 
The Higgs boson couplings also have physical boundaries in the two-dimensional parameter space of the 
2HDM (see Section 7) and hMSSM (Section 8 ) models, which are treated in a similar fashion. 

For the combination of the direct searches for invisible Higgs boson decays, confidence intervals in BRjnv 
are defined using the CF 5 procedure [45] in order to be consistent with the convention used in the indi¬ 
vidual searches. For the constraints on BRinv from the rate measurements in visible Higgs boson decay 
channels, and from the overall combination of visible and invisible decay channels, the log-likelihood 
ratio is used in order to be consistent with the convention used in deriving the Higgs boson couplings via 
the combination of visible decay channels. 

Table 1 summarises the relevant best-fit value, interval at the 68 % confidence level (CF), and/or upper 
limit at the 95% CF for physical quantities of interest. These include the overall signal strength, the scale 
factors for the Higgs boson couplings and total width, and the Higgs boson invisible decay branching 
ratio in various parameterisations. The BSM models probed with these parameters are also indicated. 
The overall signal strength measured is above 1. The extracted coupling scale factors can be similar to 
or less than 1 because the measured rate for h bb, which has a branching ratio of 57% in the SM for 
m/, = 125.36 GeV, is lower than (although still compatible with) the expected rate. 


4 Mass scaling of couplings 


The observed rates in different channels are used to determine how the Higgs boson couplings to other 
particles scale with the masses of those particles. The measurements [10] of the scale factors for the 
couplings of the Higgs boson to the Z boson, W boson, top quark, bottom quark, r lepton, and muon - 
namely [/cz, kw, Kt, Kb, Kr, x^] - are given in Model 1 of Table 1. The coupling scale factors to different 
species of fermions and vector bosons, respectively, are expressed in terms of the parameters [e, M] [46], 
where e is a mass scaling parameter and M is a “vacuum expectation value” parameter whose SM value 
is V X 246 GeV: 


Kp.i - 
KV,j = 


m2e 


(9) 
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Model 

Coupling 

Parameter 

Description 

Measurement 

j Mass scaling 

parameteri s ation 

KZ 

Z boson coupling s.f. 

[-1.06, -0.82] U [0.84,1.12] 

KW 

W boson coupling s.f. 

0.91 +0.14 

Kt 

f-quark coupling s.f. 

0.94 + 0.21 

Kb 

f7-quark coupling s.f. 

[-0.90, -0.33] U [0.28,0.96] 

Kt 

Tau lepton coupling s.f. 

[-1.22, -0.80] U [0.80,1.22] 

K^ 

Muon coupling s.f. 

< 2.28 at 95% CL 

MCHM4, 

EW singlet 

P/i 

Overall signal strength 

1 1O+0.15 

MCHM5, 

^ 2HDM Type I 

KV 

Vector boson ( W, Z) 
coupling s.f. 

1.09 + 0.07 

Kf 

Fermion (f, b,T, ...) 
coupling s.f. 

1.11 +0.16 

4 2HDM Type II, 

hMSSM 

dvM = KvIKu 

Ratio of vector boson to 
up-type fermion (f, c,...) 
coupling s.f. 

0 92+0'^ 

Kuu = K^ulKh 

Ratio of squared up-type 
fermion coupling s.f. to 
total width s.f. 

1.25 + 0.33 

^du ~ ^d!^u 

Ratio of down-type 
fermion {b,T, ...) to 
up-type fermion 
coupling s.f. 

[-1.08, -0.81] U [0.75,1.04] 

^ 2HDM 

Lepton-specific 

'^Vq = KviKq 

Ratio of vector boson to 
quark (t, b,...) 
coupling s.f. 

1 

^•^*^-0.15 

II 

Ratio of squared quark 
coupling s.f. to total 
width s.f 

^•^*^-0.20 

II 

Ratio of lepton (t, ij, e) 
to quark coupling s.f. 

[-1.34, -0.94] U [0.94,1.34] 

Higgs portal 
(Baseline config. 
of vis. & inv. 

Higgs boson 
decay channels: 
general coupling 

6 param., no 

assumption about 
Kw,z) 

KZ 

Z boson coupling s.f. 

0.99 + 0.15 

Kw 

W boson coupling s.f. 

0.92 + 0.14 

Kt 

f-quark coupling s.f. 


Kb 

/7-quark coupling s.f. 

0.61+0.28 

Kt 

Tau lepton coupling s.f. 


Kit 

Muon coupling s.f. 

< 2.25 at 95% CL 

Xg 

Gluon coupling s.f. 

0 92+0'* 

Ky 

Photon coupling s.f. 

0 90+0.16 

^Zy 

Zy coupling s.f. 

< 3.15 at 95% CL 

BRinv 

Invisible branching ratio 

< 0.23 at 95% CL 


Table 1: Measurements of the overall signal strength, scale factors (s.f.) for the Higgs boson couplings and total 
width, and the Higgs boson invisible decay branching ratio, in different coupling parameterisations, along with the 
BSM models or parameterisations they are used to probe. The measurements quoted for Models 1-5 were derived 
in Ref. [10], while those for Model 6 are derived in this paper. The production modes are taken to be the same as 
those in the SM in all cases. In Models 1-3, decay modes identical to those in the SM are taken. For Models 4—5, 
the coupling parameterisations and measurements listed do not require such an assumption, which is however made 
when deriving limits on the underlying parameters of these BSM models. No assumption about the total width is 
made for Model 6. 
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Figure 1; Two-dimensional confidence regions as a function of the mass scaling factor e and the vacuum expectation 
value parameter M. The likelihood contours where -2 In A = 2.3 and -2 In A = 6.0, corresponding approximately 
to the 68% CL (1 std. dev.) and the 95% CL (2 std. dev.) respectively, are shown for both the data and the prediction 
for a SM Higgs boson. The best ht to the data and the SM expectation are indicated as x and + respectively. 


Parameter 

Ohs. 

Exp. 

6 

M 

0.018 ±0.039 
224+J4 GeV 

0.000 ± 0.042 
246!}^ GeV 


Table 2: Observed and expected measurements of the mass scaling parameter e and the “vacuum expectation value” 
parameter M. 


where mp^i denotes the mass of each fermion species (indexed i) and myj denotes each vector-boson mass 
(indexed j). The mass scaling of the couplings, as well as the vacuum expectation value, of the SM are 
recovered with parameter values 6 = 0 and M = v, which produce Kfj = Kyj = 1. The value 6 = -1 
would correspond to light Higgs boson couplings that are independent of the particle mass. 

Combined fits to the measured rates are performed with the mass scaling factor e and the vacuum expecta¬ 
tion value parameter M as the two parameters of interest. Figure 1 shows contours of the two-dimensional 
likelihood as a function of e and M. The measured and expected values from one-dimensional likelihood 
scans are given in Table 2. The mass scaling of the couplings in the SM (e = 0) is compatible with the 
data within one std. dev. The extracted value of e is close to 0, indicating that the measured couplings to 
fermions and vector bosons are consistent with the linear and quadratic mass dependence, respectively, 
predicted in the SM. The best-fit value for M is less than v w 246 GeV because the measured overall signal 
strength /i/, is greater than 1, with the data being compatible with the SM within about 1.5 std. dev. 














5 Minimal composite Higgs model 


Minimal Composite Higgs Models (MCHM) [47-53] represent a possible explanation for the scalar nat¬ 
uralness problem, wherein the Higgs boson is a composite, pseudo-Nambu-Goldstone boson rather than 
an elementary particle. In such cases, the Higgs boson couplings to vector bosons and fermions are 
modified with respect to their SM expectations as a function of the Higgs boson compositeness scale, /. 
Corrections due to new heavy resonances such as vector-like quarks [54] are taken to be sub-dominant. 
Production or decays through loops are resolved in terms of the contributing particles in the loops, as¬ 
suming only contributions from SM particles. It is assumed that there are no new production or decay 
modes besides those in the SM. 

The MCHM4 model [47] is a minimal SO(5)/SO(4) model where the SM fermions are embedded in 
spinorial representations of SO(5). Here the ratio of the predicted coupling scale factors to their SM 
expectations, k, can be written in the particularly simple form: 

k^kv = kf= , ( 10 ) 


where ^ is a scaling parameter (with v being the SM vacuum expectation value) such that the 

SM is recovered in the limit ^ ^ 0, namely / —> cxj. The combined signal strength, jih, which is 
equivalent to the coupling scale factor, k - ^[jlh, was measured using the combination of the visible 
decay channels [10] and is listed in Model 2 of Table 1. The experimental measurements are interpreted 
in the MCHM4 scenario by rescaling the rates in different production and decay modes as functions of 
the coupling scale factors k = Ky = Kf, taking the same production and decay modes as in the SM. This 
is done in the same way as described in Section 3. The coupling scale factors are in turn expressed as 
functions of ^ using Eq. (10). 

Figure 2(a) shows the observed and expected likelihood scans of the Higgs compositeness scaling para¬ 
meter, in the MCHM4 model. This model contains a physical boundary restricting to ^ > 0, with the 
SM Higgs boson corresponding to ^ = 0. Ignoring this boundary, the scaling parameter is measured to be 
^ \ - Hh - -0.18 ± 0.14, while the expectation for the SM Higgs boson is 0 ± 0.14. The best-ht value 

observed for ^ is negative because yu/, > 1 is measured. The statistical and systematic uncertainties are of 
similar size. Accounting for the lower boundary produces an observed (expected) upper limit at the 95% 
CL of ^ < 0.12 (0.23), corresponding to a Higgs boson compositeness scale of / > 710 GeV (510 GeV). 
The observed limit is stronger than expected because ///,>! was measured [10]. 


Similarly, the MCHM5 model [48, 49] is an SO(5)/SO(4) model where the SM fermions are embedded in 
fundamental representations of SO(5). Here the measured rates are expressed in terms of ^ by rewriting 
the coupling scale factors [/cy, ] as: 


Kv = ^ 


Kf = 


l-2g 


( 11 ) 


where ^ jf^. The measurements of Ky and kf [10] are given in Model 3 of Table 1. The likelihood 

scans of ^ in MCHM5 are shown in Figure 2(b). As with the MCHM4 model, the MCHM5 model contains 
a physical boundary restricting to ^ > 0, with the SM Higgs boson corresponding to ^ = 0. Ignoring this 
boundary, the composite Higgs boson scaling parameter is determined to be ^ = -0.12 ± 0.10, while 
0.00 ± 0.10 is expected for the SM Higgs boson. As above, the best-ht value for ^ is negative because 
yU/j >1 is measured. Accounting for the boundary produces an observed (expected) upper limit at the 95% 
CL of ^ < 0.10 (0.17), corresponding to a Higgs boson compositeness scale of / > 780 GeV (600 GeV). 
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-0.8 -0.6 -0.4 -0.2 0 0.2 0.4 


-0.5-0.4-0.3-0.2 -0.1 0 0.1 0.2 0.3 


(a) MCHM4 


(b) MCHM5 


Figure 2: Observed (solid) and expected (dashed) likelihood scans of the Higgs compositeness scaling parameter, 
in the MCHM4 and MCHM5 models. The expected curves correspond to the SM Higgs boson. The line at 
-2 In A = 0 corresponds to the most likely value of ^ within the physical region ^ > 0. The line at -2 In A = 3.84 
corresponds to the one-sided upper limit at approximately the 95% CL (2 std. dev.), given ^ > 0. 


Model 

Lower limif on / 
Ohs. Exp. 

MCHM4 

710 GeV 510 GeV 

MCHM5 

780 GeV 600 GeV 


Table 3: Observed and expected lower limits at the 95% CL on the Higgs boson compositeness scale / in the 
MCHM4 and MCHM5 models. 


Figure 3 shows the two-dimensional likelihood for a measurement of the vector boson (kv) and fermion 
(kf) coupling scale factors in the [kv, kf] plane, overlaid with predictions as parametric functions of ^ for 
the MCHM4 and MCHM5 models [55-57]. Table 3 summarises the lower limits at the 95% CL on the 
Higgs boson compositeness scale in these models. 


6 Additional electroweak singlet 

A simple extension to the SM Higgs sector involves the addition of one scalar EW singlet field [41, 58- 
63] to the doublet Higgs field of fhe SM, wifh fhe doublef acquiring a non-zero vacuum expecfafion value. 
This sponfaneous symmefry breaking leads fo mixing befween fhe singlef sfafe and fhe surviving sfafe of 
fhe doublef field, resulfing in fwo CP-even Higgs bosons, where h (H) denofes fhe lighfer (heavier) of 
fhe pair. The fwo Higgs bosons, h and H, are faken fo be non-degenerafe in mass. Their couplings fo 
fermions and vecfor bosons are similar fo fhose of fhe SM Higgs boson, buf each wifh a sfrengfh reduced 
by a common scale factor, denoted by k for h and k' for H. The coupling scale facfor k (k') is fhe sine 
(cosine) of fhe h-H mixing angle, so: 

+ k'^ = \ . ( 12 ) 
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Figure 3: Two-dimensional likelihood contours in the [kv, kf] coupling scale factor plane, where -2 In A = 2.3 and 
-2 In A = 6.0 correspond approximately to the 68% CL (1 std. dev.) and the 95% CL (2 std. dev.), respectively. 
The coupling scale factors predicted in the MCHM4 and MCHM5 models are shown as parametric functions of the 
Higgs boson compositeness parameter f ! f^- The two-dimensional likelihood contours are shown for reference 

and should not be used to estimate the exclusion for the single parameter 


The lighter Higgs boson h is taken to be the observed Higgs boson. It is assumed to have the same 
production and decay modes as the SM Higgs boson does,^ with only SM particles contributing to loop- 
induced production or decay modes. In this model, its production and decay rates are modified according 
to: 

(Th = X(Th,SM 

Th = AT^XT/iSM (13) 

where o"/, denotes the production cross section, F/, denotes the total decay width, BR/,, denotes the branch¬ 
ing ratio to the different decay modes i, and SM denotes their respective values in the Standard Model. 

For the heavier Higgs boson H, new decay modes such as H ^ hh are possible if they are kinematically 
allowed. In this case, the production and decay rates of the H boson are modified wifh respecf fo fhose of 
a SM Higgs boson wifh equal mass by fhe branching rafio of all new decay modes, BR//new> as: 


o-H 

Th 

BR//,, 


k'^ xcr h,sm 


1 - BRh, new 

(l-BR//, 

new 


X F/z^sm 

) X BR//_sm,( 


(14) 


^ The decays of the heavy Higgs bosons to the light Higgs boson, for example H —> hh, are assumed to contribute negligibly to 
the light Higgs boson production rate. The contamination from heavy Higgs hoson decays (such as // —> WW) in light Higgs 
boson signal regions {h —> WW) is also taken to be negligible. 
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Upper limit on k'^ 

Obs. Exp. 

Im (U3 


Table 4: Observed and expected upper limits at the 95% CL on the squared coupling scale factor of the heavy Higgs 
boson, k’'^, in a model with an additional electroweak singlet. 


Here cthsm, r//,SM> and BR//sm,/ denote the cross section, total width, and branching ratio for a given 
decay mode (indexed i) predicted for a SM Higgs boson with mass mu- 

Consequently the overall signal strengths, namely the ratio of production and decay rates in the measured 
channels relative to the expectations for a SM Higgs boson with corresponding mass, are given by: 


(Th X BRh 
(cThX BR/,)sm 




(ThX BR// 
(cT// X BR//)sm 


(1 - BR//,new) 


(15) 


for h and H respectively, assuming the narrow-width approximation such that interference effects are 
negligible. 

Combining Eqs. (12) and (15), the squared coupling scale factor of the heavy Higgs boson can be ex¬ 
pressed in terms of the signal strength of the light Higgs boson as: 

K^ = l-Hh ■ (16) 


This equation for the squared coupling scale factor takes the same form as Eq. (10), so the same parameter 
constraints are expected. 

In particular, accounting for the lower boundary yields an observed (expected) upper limit at the 95% 
CE of k'^ < 0.12 (0.23), which is indicated in Table 4. Erom Eq. (15), this corresponds to the maximum 
signal strength for contamination by heavy Higgs boson decays in the light Higgs boson signal. Eigure 4 
shows the limits in the [////, BR//^new] plane of the heavy Higgs boson. Contours of the scale factor 
for the total width, ThI^h,sm, based on Eqs. (14) and (15), are also illustrated. These parameters are 
interesting as potential experimental observables in direct searches for heavy Higgs bosons. These results 
are independent of the mass and BR// new of the heavy Higgs boson. 


7 Two Higgs doublet model 

Another simple extension to the SM Higgs sector is the 2HDM [41, 64-66], in which the SM Higgs sector 
is extended by an additional doublet of the complex field. Eive Higgs bosons are predicted in the 2HDM: 
two neutral CP-even bosons h and H, one neutral CP-odd boson A, and two charged bosons H-. The 
most general 2HDMs predict CP-violating Higgs boson couplings as well as tree-level flavour-changing 
neufral currents. Because the latter are strongly constrained by existing data, the models considered have 
additional requirements imposed, such as the Glashow-Weinberg condition [67, 68], in order to evade 
existing experimental bounds. 
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Figure 4; Observed and expected upper limits at the 95% CL on the squared coupling scale factor, k'^, of a heavy 
Higgs boson arising through an additional EW singlet, shown in the [/UH,BR/r,new] plane. The light shaded and 
hashed regions indicate the observed and expected exclusions, respectively. Contours of the scale factor for the 
total width, F/f/rH^sM. of the heavy Higgs boson are also illustrated based on Eqs. (14) and (15). 


Both Higgs doublets acquire vacuum expectation values, vi and V 2 respectively. Their ratio is denoted 
by tan/3 = V 2 /v\, and they satisfy v\ + v\ = (246 GeV)^. The Higgs sector of the 2HDM can be 

described by six parameters: four Higgs boson masses (m/„ niH, itia, and mh^), tan/3, and the mixing 
angle a of the two neutral, CP-even Higgs states. Gauge invariance fixes the couplings of the two neutral, 
CP-even Higgs bosons to vector bosons relative to their SM values to be: 

= cos(J3-a) . 

Here W = 'W,Z and denote the SM Higgs boson couplings to vector bosons. 


The Glashow-Weinberg condition is satisfied by four types of 2HDMs [66] : 

• Type I: One Higgs doublet couples to vector bosons, while the other couples to fermions. The first 
doublet is “fermiophobic” in the limit that the two Higgs doublets do not mix. 

• Type II: This is an “MSSM-like” model, in which one Higgs doublet couples to up-type quarks 
and the other to down-type quarks and charged leptons. This model is realised in the Minimal 
Supersymmetric Standard Model (MSSM) (see Section 8). 

• Lepton-specific: The Higgs bosons have the same couplings to quarks as in the Type I model and 
to charged leptons as in Type II. 

• Flipped: The Higgs bosons have the same couplings to quarks as in the Type II model and to 
charged leptons as in Type I. 
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Coupling scale factor 

Type I 

Type II Lepfon-specific 

Flipped 

Ky 

sinfjS - a) 

Ku 

cos(a)/ sin(/3) 

Kd 

cos(a) / sinfyS) 

- sin(a)/ cosfyS) cos(a) / sin(/3) 

- sin(a)/cos(j6) 

Kt 

cos(a) / sinfyS) 

- sin(a)/cosfyS) - sin(Qr)/cosfjS) 

cos(a)/ sinfjS) 


Table 5: Couplings of the light Higgs boson h to weak vector bosons (kv), up-type quarks (/<■„), down-type quarks 
{Kd), and charged leptons (Kt), expressed as ratios to the corresponding SM predictions in 2HDMs of various types. 


Table 5 expresses the scale factors for the light Higgs boson couplings, [Ky, Ku, Kd, Ke\, in terms of a and 
tan/3 for each of the four types of 2HDMs [69]. The coupling scale factors are denoted Ky for the W and 
Z bosons, Ku for up-type quarks, Kd for down-type quarks, and Kt for charged leptons. 

The Higgs boson rate measurements in different production and decay modes are interpreted in each of 
these four types of 2HDMs, taking the observed Higgs boson to be the light CP-even neutral Higgs boson 
h. This is done by rescaling the production and decay rates as functions of the coupling scale factors 
[Ky, Ku, Kd, Ke\. The measurements of these coupling scale factors or ratios of them [10], taking the same 
production and decay modes as in the SM, are given in Models 3-5 of Table 1. These coupling scale 
factors are in turn expressed as a function of the underlying parameters, the two angles and a, using 
the relations shown in Table 5. Here the decay modes are taken to be the same as those of the SM Higgs 
boson. 

After rescaling by the couplings, the predictions agree with those obtained using the SUSHI 1.1.1 [70] 
and 2HDMC 1.5.1 [71] programs, which calculate Higgs boson production and decay rates respectively 
in two-Higgs-doublet models. The rescaled gluon fusion (ggF) rate agrees with the SUSHI prediction to 
better than a percent, and the rescaled decay rates show a similar level of agreement. The cross section 
for bbh associated production is calculated using SUSHI and included as a correction that scales with the 
square of the Yukawa coupling to the h-quark, assuming that it produces differential distributions that are 
the same as those in ggF. The correction is a small fraction of the total production rate for the regions of 
parameter space where the data would be compatible with the SM at the 95% CL. 

The two parameters of interest correspond to the quantities cos(j6 - a) and tan/3. The 2HDM possesses an 
“alignment limit” at cos(/3 - a) = 0 [66] in which all the Higgs boson couplings approach their respective 
SM values. The 2HDM also allows for limits on the magnitudes of the various couplings that are similar 
to the SM values, but with a negative relative sign of the couplings to particular types of fermions. These 
limits appear in the regions where costjS -i- a) = 0, as shown in Table 5. For example, in the Type II 
model the region where cosfjS -t a) = 0, corresponding to the sign change a -a, has a “wrong-sign 
Yukawa limit” [72, 73] with couplings similar to the SM values except for a negative coupling to down- 
type quarks. The case for the Flipped model is similar, but with a negative coupling to both the leptons 
and down-type quarks. An analogous “symmetric limit” [73] appears in the Lepton-specific model. 

Figure 5 shows fhe regions of fhe [cosfjS - Q'),fan/3] plane fhaf are excluded af a CL of af leasf 95% for 
each of fhe four types of 2HDMs, overlaid wifh fhe exclusion limifs expecfed for fhe SM Higgs secfor. 
The a and /3 paramefers are faken fo satisfy 0 < yS < njl and 0 < y6 - a < tt wifhouf loss of generalify. The 
observed and expecfed exclusion regions in cos(/3 - a) depend on fhe parficular funcfional dependence 
of fhe couplings on /3 and a, which are differenl for fhe down-fype quarks and lepfons in each of fhe 
four fypes of 2HDMs, as shown in Table 5. There is a physical boundary /cy < 1 in all four 2HDM 
fypes, fo which fhe profile likelihood ratio is resfricfed. The dafa are consisfenf wifh fhe alignmenf limif 
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at cos(/3 - a) = 0, where the light Higgs boson couplings approach the SM values, within approximately 
one std. dev. or better in each of the models. 

In each of the Type II, Lepton-specific, and Flipped models, at the upper right of the [cosfjS - a), tan/3] 
plane where tanyS is moderate, there is a narrow, curved region or “petal” of allowed parameter space with 
the surrounding region being excluded. These three allowed upper petals correspond respectively to an 
inverted sign of the coupling to down-type fermions (tau lepton and bottom quark), leptons (t and //), or 
the bottom quark. These couplings are measured with insufficient precision to be excluded. There is no 
upper petal at high tan/3 in Type I as all the Yukawa couplings are identical. 

In each of the four 2HDM types a similar petal is possible at the lower right of the [cosfyS- a), tan /S] plane. 
For the Type I, Type II, Lepton-specific, and Flipped models, this lower petal corresponds respectively to 
an inverted coupling to fermions, up-type quarks, all quarks, and lastly the up-type quarks and leptons. 
In all four cases, the lower petal is rejected since an inverted top quark coupling sign is disfavoured. The 
top quark coupling is extracted primarily through its dominant effect in ggF Higgs production, as well as 
by resolving the Higgs boson decays to diphotons, with one contribution being from the top quark. 

For this analysis, only the range 0.1< tan^S <10 was considered. The regions of compatibility extend to 
larger and smaller tan/3 values, but with a correspondingly narrower range of cos(/3 - a). The confidence 
intervals drawn are derived from a distribution with two parameters of interest, corresponding to 
the quantities cosfjS - a) and tan/3. However, at cosfyS - a) = 0 the likelihood is independent of the 
model parameter /3, effectively reducing the number of parameters of interest locally to one. Hence 
the test-statistic distribution for two parameters of interest that is used leads to some overcoverage near 
cos(/3 - cr) = 0. 


8 Simplified Minimal Supersymmetric Standard Model 


Supersymmetry provides a means to solve the hierarchy problem by introducing superpartners of the 
corresponding SM particles. Many supersymmetric models also provide a candidate for a dark-matter 
particle. In the Minimal Supersymmetric Standard Model [74-80], the mass matrix of the neutral CP- 
even Higgs bosons h and H can be written as [81]: 


m| cos^/3 - 1 - sin^/3 

-(m| - 1 - m^) sin/3 cos p 

_ 1 _ 


AM^^' 

-(m| + m^) sin/3 cos/3 

/n| sin^yS -I- cos^/3 



AAI22 


with radiative corrections being included through the 2x2 matrix AAl^^. 

A simplified approach to the study of the MSSM Higgs sector, known as the hMSSM [81-83], consists 
of neglecting the terms AAIjj and AAIj 2 - The remaining term AAI 22 , which contains the dominant 
corrections from loops involving top quarks and stop squarks, is traded for the lightest mass eigenvalue 
m/,. The scale factors for the Higgs boson couplings to vector bosons, up-type fermions, and down-type 
fermions ([/cy, /<■„, Kd\), can be expressed as functions of the free parameters [wta, tan/3] (in addition to nih) 
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(b) Type II 
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Figure 5; Regions of the [cos(/3 - a), tan p] plane of four types of 2HDMs excluded by fits to the measured rates of 
Higgs boson production and decays. The likelihood contours where -2 In A = 6.0, corresponding approximately to 
the 95% CL (2 std. dev.), are indicated for both the data and the expectation for the SM Higgs sector. The cross in 
each plot marks the observed best-fit value. The light shaded and hashed regions indicate the observed and expected 
exclusions, respectively. The a and p parameters are taken to satisfy 0 < p < njl and 0 < p - a < n without loss 
of generality. 
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as [81-83]: 


_ i’^;(m/i,tan/l)+tany8 i„(OTA,tanj6) 

1 +tan2 p 

-J 1 +tan2 p 
Ku = 5i,(mA, tanyS)—— 

Kd = SdiniA, tan/3) ^1 +tan^/3 , 


where the functions and Sd are given by: 


Su = 


Sd = 


, 1 +--- 

(m^ + m|) tan/3 
tan^/i-m^(l+tan^/3) 


( 18 ) 


(19) 


and mz is the mass of the Z boson. 

To test the hMSSM model, the measured production and decay rates are expressed in terms of the cor¬ 
responding coupling scale factors for vector bosons (/ty), up-type fermions (a:„), and down-type fermions 
{Kd)- The observed Higgs boson is taken to be the light CP-even neutral Higgs boson h. In the hMSSM, 
it is assumed to have the same production and decay modes as in the SM. For comparison. Model 4 of 
Table 1 lists the measurements of ratios of the coupling scale factors [/cy, Ku, Kd) [10]. The coupling scale 
factors are in turn cast in terms of tua and tanyS using Eq. (18). A correction is applied for bbh associated 
production as a function of the /7-quark Yukawa coupling as described in Section 7. 

Loop corrections from stops in ggF production, which can decrease the rate by 10-15% for a light 
stop [84], and in diphoton decays are neglected. Light tau sleptons (staus) with large mixing could 
enhance the diphoton rate by up to 30% at tan/3 = 50 [84], and charginos could modify the diphoton rate 
by up to 20% [74, 75, 77, 85]; these effects are not included in the hMSSM model. 

Additional corrections in the MSSM would break the universality of down-type fermion couplings, res¬ 
ulting for example in Kb + Kt These are generally sub-dominant effects [81-83] and are not included. 
The MSSM includes other possibilities such as Higgs boson decays to supersymmetric particles, decays 
of heavy Higgs bosons to lighter ones [86], and effects from light supersymmetric particles [84], which 
are not investigated here. This model is therefore not fully general but serves as a useful benchmark, 
particularly if no direct observation of supersymmetry is made. 

Contours of the two-dimensional likelihood in the [wta, tanyS] plane for the hMSSM model are shown in 
Figure 6. The data are consistent with the SM decoupling limit at large tua- The observed (expected) 
lower limit at the 95% CL on the CP-odd Higgs boson mass is at least niA > 370 GeV (310 GeV) for 
1 < tanyS < 50, increasing to 440 GeV (330 GeV) at tan/3 = 1. The observed limit is stronger than 
expected because the measured rates in the h ^ yy{2V\ (expected to be dominated by a IT boson loop in 
the SM) and h ZZ* —> [22] channels are higher than predicted by the SM, but the hMSSM model 
has a physical boundary /ry < 1 so the vector-boson coupling cannot be larger than the SM value. The 
physical boundary is accounted for by computing the profile likelihood ratio with respect to the maximum 
likelihood obtained within the physical region of the parameter space, mA > 0 and tan/3 > 0. The region 
1 < tanyS < 50 is shown; at significantly smaller or larger values of tan/3, the hMSSM model is not a 
good approximation of the MSSM. For tan/3 < 1, the couplings to SM particles receive potentially large 
corrections related to the top sector that have not been included [87]. 


17 













[GeV] 


Figure 6; Regions of the [niA, tan^S] plane excluded in the hMSSM model via direct searches for heavy Higgs 
bosons and hts to the measured rates of observed Higgs boson production and decays. The likelihood contours 
where -2 In A = 6.0, corresponding approximately to the 95% CL (2 std. dev.), are indicated for the data (solid 
lines) and the expectation for the SM Higgs sector (dashed lines). The light shaded or hashed regions indicate the 
observed exclusions. The SM decoupling limit is ^ oa. 


The constraints in the [ma, tan/3] plane of the hMSSM model from various direct searches for heavy 
Higgs bosons are also overlaid in Figure 6. The constraints from the following searches are shown. 

• H/A TT search via both ggF and bbh associated production [88]. 

• Heavy CP-odd Higgs boson A produced via ggF and decaying to Zh with Z ^ ee, jufi, or vv and 
h bb [89]. 

• Heavy CP-even Higgs boson H produced via ggF and decaying to WW —> bv£v, (vqq [90] or 
ZZ ^ 4i, iiqq, i£bb, ttw [91] (I - e, yc, q - u, d, c, s). 

• Charged Higgs boson H- production in association with a top quark [92]. 

The cross sections for ggF and bbh associated production in the five-flavour scheme hMSSM model have 
been calculated with SUSHI 1.5 . Q [70]. The calculation for ggF includes the complete massive top and 
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bottom loop corrections at next-to-leading-order (NLO) QCD [93], the top quark loop corrections in the 
heavy-quark limit of QCD at next-to-next-to-leading-order (NNLO) [94-96], and EW loop corrections 
due to light quarks up to NLO [97, 98]. The bbh associated production in the five-flavour scheme includes 
corrections up to NNLO in QCD [99]. The production of heavy Higgs bosons has also been calculated 
in the four-flavour scheme at NLO in QCD [100, 101], and the result has been combined with the five- 
flavour scheme using an empirical mafching procedure [102]. The branching rafios have been calculated 
using HDECAY 6.4.2 [103]. 


9 Probe of invisible Higgs boson decays 

9.1 Direct searches for invisible decays 

Linal sfafes wifh large missing fransverse momenfum associafed wifh lepfons or jefs offer fhe possibilify 
of direcf searches for h invisible [104-11 1]. In fhese searches, no excess of evenfs was found and upper 
limifs were sef on fhe Higgs boson producfion cross secfion fimes fhe branching rafio for h —> invisible 
decays. Assuming fhaf fhe Higgs boson producfion cross sections and accepfances are unchanged relafive 
fo fhe SM expecfafions, upper bounds on fhe branching rafio of invisible Higgs boson decays, BRinv, 
were obfained from fhe cr x BR measuremenfs. The ATLAS and CMS Collaborafions sef upper limifs af 
fhe 95% CL of 28% [32] and 65% [112], respectively, on fhe branching rafio for invisible Higgs boson 
decays by searching for vecfor-boson fusion production of a Higgs boson fhaf decays invisibly. Using fhe 
Zh ^ signafure, weaker bounds were obfained by bofh ATLAS and CMS, giving upper limifs 

of 75% [33] and 83% [112], respectively. By combining searches in Z{€l)h and Z(bb)h, CMS obfained 
an upper limif of 81% [1 12]. A combination of fhe searches in VBL and Zh producfion was carried ouf by 
CMS, giving a combined upper limif of 58% [112]. Using fhe associated producfion wifh a vector boson, 
Vh, where V - W or: Z, V ^ jj, and h —> invisible, ATLAS sef an upper bound of 78% [34]. Ofher 
searches for invisible Higgs boson decays in evenfs wifh large in associafion wifh one or more jefs 
were also performed [113-116], buf fhese searches are less sensitive fo Higgs-mediafed interactions. In 
fhe SM, fhe process h —> ZZ* —> 4y is an invisible decay mode of fhe Higgs boson, buf fhe branching 
rafio is 0.1% [41], which is below fhe sensitivities of fhe aforementioned direcf searches. 

A sfafisfical combinafion of fhe following direcf searches for invisible Higgs boson decays is performed: 

(1) The Higgs boson is produced in fhe VBL process and decays invisibly [32]. The signafure of fhis 
process is fwo jefs wifh a large separation in pseudorapidify, forming a large invarianf dijef mass, 
fogelher wifh large 

(2) The Higgs boson is produced in associafion wifh a Z boson, where Z ^ tl and fhe Higgs boson 
decays fo invisible particles [33]. The signafure in fhis search is fwo opposite-sign and same-flavour 
lepfons (elecfrons or muons) wifh large missing fransverse momenfum. 

(3) The Higgs boson is produced in associafion wifh a vector boson V (IT or Z), where V ^ H and 
fhe Higgs boson decays fo invisible particles [34]. The signafure in fhis search is fwo jefs whose 
invarianf mass m^j is consisfenf wifh fhe V mass, fogefher wifh large missing fransverse momenfum. 

To combine fhe measuremenfs, fhe searches need fo be performed in non-overlapping regions of phase 
space or fhe combinafion musf accounf for fhe overlap in phase space. The Zh, search does nol 

overlap wifh fhe ofher searches for h invisible because a vefo on evenfs confaining jefs was required. 
The overlap due fo possible inefficiency in fhe vefo requiremenfs is negligible. The VBL —> jj + L™®® 
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and the WjZh jj + searches also do not overlap in their phase spaces because the former requires 
a large dijet invariant mass (above rujj > 500 GeV) and latter imposes the requirement that the dijet 
invariant mass must be consistent with the associated vector boson mass within 50 < nijj < 100 GeV 
and imposes a veto on forward jets. The same overlap removal requirements were applied in data to both 
the signal and control regions in the various searches, making the control regions used for background 
estimation non-overlapping. 

The following nuisance parameters are treated as being fully correlated across the individual searches, 
with the rest being uncorrelated: 

• Uncertainty in the luminosity measurements. This impacts the predicted rates of the signals and the 
backgrounds that are estimated using Monte Carlo simulation, namely ggF, VBF, and Vh signals, 
and tt, single top, and diboson backgrounds. 

• Uncertainties in the absolute scale of the jet energy calibration and on the resolution of the jet 
energy calibration. 

• Uncertainties in the modelling of the parton shower. 

• Uncertainties in the renormalisation and factorisation scales, as well as the parton distribution func¬ 
tions. This affects the expected numbers of signal events in the ggF, VBF and Vh production chan¬ 
nels. 

The uncertainty in the soft component of the missing transverse momentum has a significant impact in 
the VF/Z/i ^ jj + F™®* channel. Its impact is much smaller in the other searches and not included as a 
nuisance parameter. This uncertainty is therefore not correlated across all the searches. 

The limit on the branching ratio oih ^ invisible, defined in Eq. (3), is compufed assuming fhe SM pro¬ 
duction cross secfions of fhe Higgs boson. This is done using a maximum-likelihood fif fo fhe evenf counfs 
in fhe signal regions and fhe dafa confrol samples following fhe CL 5 modified frequenfisf formalism wifh 
a profile likelihood-rafio fesf stafisfic [40]. 

The sfafisfical combination of direcf searches for invisible Higgs boson decays resulfs in an observed 
(expected) upper limif af fhe 95% CL on BRmv of 0.25 (0.27). Figure 7 shows fhe scan of fhe CL as a 
function of BR(/j ^ invisible) for fhe sfafisfical combination of direcf searches. The limif obfained wifh 
fhe CL 5 mefhod is consisfenf fo fwo significanf figures wifh fhe limif based on fhe log likelihood rafio. 
Table 6 summarises fhe limifs from direcf searches for invisible Higgs boson decays and fheir sfafisfical 
combination. The combined limif is dominafed by fhe VBL ^ jj + channel, which is by far fhe 
mosf sensitive. 


9.2 Combination of visible and invisible decay channels 

The measurements of Higgs boson visible decay rates are complementary to the direct searches for invis¬ 
ible decays since they are indirectly sensitive to undetectable decays as well. The visible decay rates are 
used to extract the sum of the branching fractions to invisible and undetectable final sfafes. A conservafive 
limif on fhe invisible branching rafio is fhen inferred by assuming fhe undefecfable branching rafio fo be 
negligibly differenl from fhe SM expecfafion of approximately zero. Lor example, a significanf excess in 
a search for Higgs boson decays fo a new particle would further tighten the indirect limit on the invisible 
branching ratio. It would not affect the limits from direct searches for invisible Higgs boson decays. 
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Figure 7; The (1 - CL) versus BR(/! —> invisible) scan for the combined search for invisible Higgs boson decays. 
The horizontal dashed lines refer to the 68% and 95% confidence levels. The vertical dashed lines indicate the 
observed and expected upper bounds at the 95% CL on BR(/i —»invisible) for the combined search. 


Channels 

Ohs. 

Upper limit on BR(/z 
-2 std. dev. -1 std. dev. 

^ inv.) at the 95% CL 
Exp. -1-1 std. dev. 

+2 std. dev. 

VBF/i 

0.28 

0.17 

0.23 

0.31 

0.44 

0.60 

z(^ a)h 

0.75 

0.33 

0.45 

0.62 

0.86 

1.19 

V(^ jj)h 

0.78 

0.46 

0.62 

0.86 

1.19 

1.60 

Combined Results 

0.25 

0.14 

0.19 

0.27 

0.37 

0.50 


Table 6: Summary of upper bounds on BR(/i —> invisible) at the 95% CL from the individual searches and their 
combination. The Higgs boson production rates via VBF and Vh associated production are assumed to be equal to 
their SM values. The numerical bounds larger than 1 can be interpreted as an upper bound on ct/ctsm, where ctsm 
is the Higgs boson production cross section in the SM. 


The overall upper limit on the branching ratio of the Higgs boson to invisible final states, BRinv, is derived 
using a statistical combination of measurements from both the visible and invisible Higgs boson decays. 
The visible decay channels are h^yy, h^ZZ*^4£, WW* ^ €v£v, h^Zy, h ^ tt, h ^ jufi, and 

h bb, with a variety of production mode selections used. The invisible decay channels are described in 
Section 9.1 and involve the Higgs boson being produced via VBF or Z{££)h, and then decaying invisibly. 
The V{jj)h production mode is not included due to overlap of the event selection with the 0-lepton 
category of the Vh{bb) measurement. 

The extraction of BRjnv is performed using a coupling parameterisation that includes separate scale factors 
for the couplings of the Higgs boson to the W boson, Z boson, top quark, bottom quark, tau lepton, and 
muon, as well as scale factors for effective loop-induced couplings to gluons, photons, and Zy to absorb 
the possible contributions of new particles through loops. The Higgs boson production modes are taken 
to be the same as those in the SM. 

As for the visible decay rates alone, the invisible branching ratio is conservatively estimated by taking the 
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undetectable branching ratio to be zero. Thus the scale factor is equal to the ratio of the total width of 
the Higgs boson to the SM expectation, T/j/T/j sm> and can be expressed in terms of BRinv as: 

kI = Th/Th,sM = Z Af^BRy/Cl - BRinv) ■ ^20) 

The production and decay rates of all channels are fit with functions of [kw, kz, Kt, Kb, Kt, Kf^, Kg, Ky, Kzy, 
BRinv]. Each parameter of interest is determined by treating the others as nuisance parameters. 

With the visible decay channels alone, the assumption that the vector-boson coupling is less than or 
equal to the SM expectation {kv < 1) produced an observed (expected) upper limit at the 95% CL of 
BRinv < 0.49 (0.48) [10]. With only the invisible decay channels (including V{jj)h), the assumption that 
the vector-boson and gluon couplings are identical to their SM expectations (/c, = 1) yields an upper limit 
at the 95% CL of 0.25 (0.27), as described in Section 9.1. 

When the visible and invisible decay channels (but excluding V{jj)h) are combined, no assumption about 
Ky or Other couplings is made beyond that on the undetectable branching ratio. The resulting observed 
likelihood scan as a function of BRinv is shown in Ligure 8. The fitted values of the Higgs boson couplings 
are similar to those in the SM, so the measured overall signal strength(/r/; >1) is accommodated in the fit 
by a best-fit value of BRinv that would be negative. This decreases the Higgs boson total width, which is 
inversely proportional to the signal strength. Scans of the observed likelihoods of the invisible and visible 
decay channels separately are also shown. 

The observed (expected) upper limit at the 95% CL is BRinv < 0.23 (0.24) using the combination of all 
channels and accounting for the physical boundary.^ This baseline configuration, with a maximally gen¬ 
eral set of independent coupling parameters and no explicit assumption about the value of xvyz. provides 
the main result for the combination of invisible and visible decay channels. The results are dominated by 
the direct searches for invisible decays. The expected limit demonstrates a fractional improvement of 11 % 
in sensitivity compared to the invisible decay channels alone. In addition, it is more model-independent 
than the combination of invisible (visible) decay channels alone, because it does not assume the vector- 
boson couplings to be equal to (less than or equal to) their SM values. The absolute couplings to the Z 
boson, W boson, top quark, bottom quark, tau lepton, muon, gluon, photon, and Zy are extracted in the 
baseline scenario. The measurements or limits are given, along with the upper limit on BRinv, in Model 6 
of Table 1 . 

As an alternate scenario, if the assumption /cvyz < 1 is added for the combination of channels, the observed 
(expected) limit is 0.23 (0.23). This allows a useful comparison with the results of the invisible or visible 
decay channels alone, which apply a similar assumption. The results from the invisible and visible decay 
channels separately, as well as their statistical combination, are summarised in Table 7 for the coupling 
parameterisation and assumptions about xvyz used. The baseline results for the combination of invisible 
and visible decays are highlighted. 

A less general coupling parameterisation [/cy, Kf, Kg, Ky, Kzy, BRinv] was also considered. The reduction 
in the number of degrees of freedom would make the limit signihcantly more model-dependent and only 
improve the sensitivity marginally, so this parametrisation was not used. 


The observed upper limit at the 90% CL is BRinv < 0.22, which is used in Section 9.3. 
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Figure 8: Observed likelihood scans of the Higgs boson invisible decay branching ratio using direct searches for 
invisible Higgs boson decays, rate measurements in visible Higgs boson decays, and the overall combination of 
invisible and visible decay channels. The line at -2 In A = 0 corresponds to the most likely value of BRinv within 
the physical region BRinv S 0. The line at -2 In A = 3.84 corresponds to the one-sided upper limit at approximately 
the 95% CL (2 std. dev.), given BRinv ^ 0. 


Decay channels 

Coupling parameterisation 

Ki assumption 

Upper limit on BRinv 
Obs. Exp. 

Invisible decays 

Kz-, A'f) ^Zyf ^^inv] 

KW,Z,g - 1 

0.25 

0.27 

Visible decays 

KZi ^Zyi ^^^inv] 

KW,Z ^ 1 

0.49 

0.48 

Inv. & vis. decays 

®®^inv] 

None 

0.23 

0.24 

Inv. & vis. decays 

[a'^^, KZi ^bf ^g ^yi ^Zy^ ^^^inv] 

KW,Z ^ 1 

0.23 

0.23 


Table 7: Summary of upper limits on BR(/r —» invisible) at the 95% CL from the combination of direct searches 
for invisible Higgs boson decays, the combination of measurements of visible Higgs boson decays, and the overall 
combination using both the invisible and visible Higgs boson decays. The results are derived using different as¬ 
sumptions about Kw,z- The results with the baseline conhguration for the combination of invisible and visible decay 
channels are indicated in bold. 
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9.3 Higgs portal to dark matter 


Many “Higgs portal” models [18, 19, 33, 62, 104, 113, 117, 118] introduce an additional weakly interact¬ 
ing massive particle (WIMP) as a dark-matter candidate. It is taken to interact very weakly with the SM 
particles, except for the Higgs boson. In this study, the coupling of the Higgs boson to the WIMP is taken 
to be a free parameter. The Higgs portal model, where the Higgs boson is taken to be the only mediator, 
is a special case of the spin-independent limits obtained by direct detection experiments. 

To compare with these direct searches, the observed upper limit combining all the visible and invisible 
Higgs boson decay channels using the most general baseline parameterisation is calculated at the 90% 
CL: this gives BRinv < 0.22. This limit is translated into constraints on the coupling of the WIMP to the 
Higgs boson as a function of its mass [19]. This is done for WIMP masses less than half the Higgs boson 
mass, under the assumption that the resulting Higgs boson decays to WIMP pairs account entirely for 
BRinv. Any additional contributions to BRinv from other new phenomena would produce a more stringent 
limit. The partial width for Higgs boson decays to a pair of dark-matter particles depends on the spin of 
the dark-matter particle. It is given for scalar, Majorana fermion, or vector dark-matter candidates (where 
the Majorana fermion is motivated by neutralinos in supersymmetry) as: 


scalar 5 : r"''(/j ^SS) = 

i2 

'^hSS 

fermion / : F‘"''(/j —> //) = 

^Iff 

A2 

vector V : r'"'(/j ^ VV) - 

i2 

'^hVV 




,2o3. 


647r 


Sllnniy 


mr 


X 1 -4^+12- 


m. 


mr 


ml 


( 21 ) 


Here Ahss^ and Ahvv are the couplings of the Higgs boson to dark-matter particles of corres¬ 

ponding spin, V denotes the vacuum expectation value of the Higgs boson, and = .yjl - Am^ltri^ is a 
kinematic factor associated with the two-body h —> XX decay (x = S, P, or /). For the cases of a fermion 
or vector boson WIMP in this effective field theory approach, the new physics scale A is assumed to be 
at the TeV scale or higher, well above the probed scale at the SM Higgs boson mass. These equations are 
used to deduce the coupling of the Higgs boson to the WIMP for each of the three possible WIMP spins. 
The coupling is then re-parameterised in terms of the cross section for scattering between the WIMP and 
nucleons via Higgs boson exchange, cr^-v? [19]: 


scalar S : 


fermion / : 


vector V : 


4 r2 

= r2 ^NfN 

^ 167rm^(m5-I-mA?)2 

_ 4ff KfWf 

^ ^ Anm^{mf + 

4 

- J2 ^NfN 

^ hvv _|. rriN)^ 


( 22 ) 


where ~ 0.94 GeV is the nucleon mass, and = 0.33];QQy is the form factor associated with the 
Higgs boson-nucleon coupling, computed using lattice QCD [19]. Upper limits at the 90% CL on the 
WIMP-nucleon scattering cross section are derived as a function of the WIMP mass as shown 
in Figure 9. The hashed bands indicate the uncertainty resulting from the systematic variation of the form 
factor fy. They are compared with limits from direct searches for dark matter [119-127] at the confidence 
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levels indicated. The ATLAS limits on the WIMP-nucleon scattering cross section are proportional to 
those on the invisible decay branching ratio, as evident from Eqns. (21)-(22). They are weaker (stronger) 
at low mass for scalar (Majorana and vector) WIMPs, and degrade as approaches ni/,/2 as expected 
from kinematics. The limits are shown for > 1 GeV, but extend to WIMP masses smaller than this 
value. 


10 Conclusions 

Higgs boson coupling measurements from the combination of multiple production and decay channels 
(/j—> 77 , ZZ*—> 4^, WW*^ €viv, Zy, h bb, h ^ tt, h ^ fiju) have been used to indirectly 
search for new physics. The results are based on up to 4.7 fb“^ of pp collision data ai ^fs - 1 TeV and 
20.3 fb“^ af ^[s - % TeV recorded by the ATLAS experiment at the LHC. No significant deviation from 
the SM expectation is found in the observables studied, which are used to constrain various models of 
new phenomena. 

The mass dependence of the couplings is consistent with the predictions for a SM Higgs boson. Con¬ 
straints are set on Minimal Composite Higgs Models, models with an additional electroweak singlet, and 
two-Higgs-doublet models. A lower limit at the 95% CL is set on the mass of a pseudoscalar Higgs boson 
in the hMSSM of 370 GeV. Results from direct searches for heavy Higgs bosons are also interpreted in 
the hMSSM. In addition, direct searches for invisible Higgs boson decays in the VBF, Z{£()h, and V{jj)h 
production modes are combined to set an upper bound at the 95% CL on the Higgs boson invisible de¬ 
cay branching ratio of 0.25. Including the coupling measurements in visible decays further improves the 
upper limit to 0.23. The limit on the invisible decay branching ratio is used to constrain the rate of dark 
matter-nucleon scattering in a model with a Higgs portal to dark matter. 
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Figure 9 ; ATLAS upper limit at the 90 % CL on the WIMP-nucleon scattering cross section in a Higgs portal model 
as a function of the mass of the dark-matter particle, shown separately for a scalar, Majorana fermion, or vector- 
boson WIMP. It is determined using the limit at the 90 % CL of BRinv < 0.22 derived using both the visible and 
invisible Higgs boson decay channels. The hashed bands indicate the uncertainty resulting from varying the form 
factor fff by its uncertainty. Excluded and allowed regions from direct detection experiments at the confidence 
levels indicated are also shown [ 119 - 127 ]. These are spin-independent results obtained directly from searches for 
nuclei recoils from elastic scattering of WIMPs, rather than being inferred indirectly through Higgs boson exchange 
in the Higgs portal model. 
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